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Abstract 

This paper outlines our recent attempts to model the growth and form of micro- 
bialites from the perspective of the statistical physics of evolving surfaces. Micro- 
bialites arise from the environmental interactions of microbial communities (mi- 
crobial mats). The mats evolve over time to form internally laminated organosed- 
imentary structures (stromatolites). Modern day stromatolites exist in only a few 
locations, whereas ancient stromatolitic microbialites were the only form of life for 
much of the Earth's history. They existed in a wide variety of growth forms, ranging 
from almost perfect cones to branched columnar structures. The coniform structures 
are central to the heated debate on the oldest evidence of life. We proposed a biotic 
model which considers the relationship between upward growth of a phototropic or 
phototactic biofilm and mineral accretion normal to the surface. These processes 
are sufficient to account for the growth and form of many ancient stromatolities. 
These include domical stromatolites and coniform structures with thickened apical 
zones typical of Conophyton. More angular coniform structures, similar to the stro- 
matolites claimed as the oldest macroscopic evidence of life, form when the photic 
effects dominate over mineral accretion. 



1 Introduction 



The Kardar, Parisi and Zhang (KPZ) equation [1] describing the temporal and 
spatial evolution of an interface has found widespread applicability in a diverse 
range of phenomena [2]. These include the description of evolving fronts in ex- 
periments on paper burning, chemical electrodeposition and yeast colonies. 
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KPZ originally suggested similarities with geological stratification. Now stro- 
matolites are internally laminated biosedimentary structures [3,4] produced as 
a consequence of some of the environmental interactions of benthic microbial 
communities (BMC) [5]. These microbiahtes [6] have been described as Earth's 
default ecosystem [7] , re-establishing themselves when other life- forms are ex- 
tinguished by environmental stress [8]. Stromatolitic microbiahtes provide an 
almost continuous record of the subsequent evolution of environmental condi- 
tions on Earth. They flourish wherever and whenever conditions favour their 
mineralisation, and where they are neither out-competed by faster-growing 
higher organisms nor consumed by grazers. 

Most importantly, stromatolitic microbiahtes preserve the only macroscopic 
evidence of life prior to the appearance of macro-algae [9] . The biogenicity of 
stromatohtes older than 3,200 miUion years is unclear [10,11,12,13]. If they 
are indeed biotic, they are the oldest morphological evidence for hfe, now 
that the identification of 3,300 to 3,500 million year old microfossils [14] has 
been challenged [15,16]. Central to this debate are the sharply peaked conical 
stromatolites described from the 3,450 million year old Warrawoona Group in 
Western Austraha [13]. 

We recently proposed a KPZ-like model for stromatolite morphogenesis that 
endorses a biotic origin for coniform stromatolites [17]. The model describes 
interaction between upward growth of a phototropic or phototactic microbial 
mat and mineral accretion normal to the surface of the mat. Domical struc- 
tures are formed when mineral accretion dominates. When vertical growth 
dominates, coniform structures evolve that reproduce the features of Cono- 
phyton, a stromatolite that flourished in certain low-sedimentation environ- 
ments for much of the Proterozoic [18]. Increasing the dominance of vertical 
growth produces sharply-peaked structures similar to those described in the 
Warrawoona Group. Here we consider this model in (2-|-l)-dimensions and 
give some preliminary results on the solutions to the model equation. 



2 The model 

The proposed model [17] for stromatolite morphogenesis involves two processes 
only: (i) upward growth of a phototropic or phototactic BMC, (ii) mineral 
accretion normal to the surface. Here the function h{x, y, t) represents the 
height of the profile above a horizontal baseplane which evolves in time t 
according to the equation 
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We interpret v as the average rate of vertical growth due to photic response of 
microbes and A as the average rate of surface-normal growth due to mineral 
accretion. 

If only radially symmetric solutions are considered then the problem reduces 
to the (l-l-l)-dimensional case after transforming to polar coordinates. Our 
approach to use the KPZ equation to describe stromatolites [17,19,20,21] dif- 
fers from others [22,23] in that we interpret our model in purely biological 
terms. We also disregard the usual KPZ noise term, thus allowing us to write 
down explict solutions to the model equation. 



3 Preliminary results and discussion 



Although non-linear, equation (1) can be solved with a change of variables 
using the method of characteristics [19] and prescribed initial profiles. The 
choice of initial profile is important. It is worth noting that cone-like profiles 
arise naturally in deformations of thin flat sheets [24] . For simphcity, we con- 
sider here a radially symmetric cone centred on the origin with base radius R 
and initial height /i(0, 0,0) — ho-^ Then the full solution to equation (1) is 
given by 



h{x,y,t) 



^0 - ^ + 



y + X)t ifr<^, 

if r>i?+^, 



{v + \)t 



(2) 



where r = \/a?~+~y^- From the solution we see that the convex sharp protru- 
sion at the tip of the initial cone gives rise to a smooth radially symmetric 
paraboloid centred on the origin, expanding outwards with time. The bound- 
ary of the cone's base represents a concave discontinuity in the spatial deriva- 
tive, and gives rise to a shock which propagates outwards with time, effectively 
expanding the base of the cone. 

In a similar way one can consider more complex cases, such as a paraboloid 
with an elliptical base. The paraboloid can be either convex or concave, the 
degree to which can be tuned by an adjustable parameter a. However, in 
all the cases we have studied, an isolated protrusion centred at the origin 
and evolving under equation (1) approaches the smooth radially symmetric 
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paraboloid solution, 

2 

h{x, y, t) = /i(0, 0> 0) - ^ + + >^)t- (3) 

For an initial profile consisting of a field of protrusions, the long-term solution 
will therefore consist of smooth radially symmetric paraboloids separated by 
shocks. Indeed, a feature of the KPZ equation is that the long-term solutions 
are radially symmetric, independent of the initial profiles. 

The parameter A controls the rate at which the shape of the paraboloid 
changes, as well as the rate of vertical growth (via the {v + X)t term). The 
parameter v, on the other hand, affects only the rate of vertical growth. The 
relative magnitudes of A and v therefore determine the rate at which the 
paraboloid changes in shape relative to its rate of upward growth. For exam- 
ple, if A is very small compared to v, then even after a long period of time, 
as measured by the amount of upward growth, the paraboloid will still closely 
resemble its initial form. If A is comparable to or greater than v, the tran- 
sition from the initial profile to a smooth radially symmetric paraboloid will 
take place over a much shorter period of vertical growth. These regimes are 
illustrated in the cross-sections shown in Figure 1, which is to be compared 
with Figure 1 of Ref. [17] obtained for the (l-l-l)-dimensional case. 

The results provide possible explanations for variations in coniform stroma- 
tolite morphogenesis. Together with field evidence, they support the inter- 
pretation that the vertical growth parameter v represents photic response of 
the BMC rather than sediment deposition. If the converse were true, coniform 
stromatolites would only form under conditions of high sedimentation which is 
precisely contrary to field evidence. Indeed, while sediment deposition would 
tend towards the smoothing of surface irregularities, growth due to photic re- 
sponse would tend to accentuate them. Our model shows that a combination of 
vertical phototropic or phototactic microbial growth and surface-normal min- 
eral accretion can produce coniform forms and structures analogous to those 
found in both Archaean and Proterozoic coniform stromatolites. For exam- 
ple, there is a striking similarity between the model forms shown in Figure 1 
and the sharply-peaked coniform stromatolites in the Warrawoona Group [13], 
thus supporting their biogenic origin and reinforcing the probability that pho- 
tosynthetic microbes were components of Archaean BMCs. 

The various cases modelled in Figure 1 can all be matched in Proterozoic 
Conophytons. Figure 2 shows a comparison between Conophytons observed 

in the Northern Territory of Australia and a longitudinal cross-section of a 
preliminary simulation of two evolving contiguous cones. Figure 3 shows a 
comparison between Conophytons and a transverse cross-section of an evolving 
paraboloid field. Our aim here is to determine the environmental parameters 
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responsible for the various Conophyton growth forms. 

We beheve that our results also shed some light on why, after flourishing for 
much of the Proterozoic, Conophytons virtually disappeared in the Neopro- 
terozoic [25]. This demise has been linked to evolutionary changes in BMCs 
[23], but since these would not have limited photic response, this seems un- 
tenable. Conophytons represent an effective growth strategy that is especially 
vulnerable to predation and competition [26] . Perhaps their demise is best ex- 
plained in terms of the evolution of greater biological diversity in their quiet 
marine environments. 
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Fig. 1. Equal time snapshots of the surface height profile for different values of the 
initial shape parameter a and the growth parameters A and v. First row: a = 20, 
second row: a = 1.5, third row: a = —1. First column: A = 0.01, v = 1.99, second 
column: A = 0.1, f = 1.9, third column: A = l,f = 1. 
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Fig. 2. (left) 1,700 million year old Conophytons from the Dungaminnie Forma- 
tion, Northern Territory, Australia. Coin 28mm for scale, (right) Longitudinal 
cross-section of a preliminary simulation that does not attempt to take the complex 
paraboloid into account. 




Fig. 3. (left) 1,670 million year old Conophytons from the Macarthur Group, North- 
ern Territory, Australia. Lens cap 25mm for scale, (right) Tansverse cross-section of 
a preliminary simulation. 



8 



